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Chapter 1 

Wave turbulence in 
magnet ohydro dynamics (MHD) 



This chapter reviews the recent progress made mainly during the last two decades 
on wave turbulence in magnetized plasmas (MHD, Hall MHD and electron MHD) in 
the incompressible and compressible cases. The emphasis is made on homogeneous 
and anisotropic turbulence which usually provides the best theoretical framework 
to investigate space and laboratory plasmas. The interplanetary medium and the 
solar atmosphere are presented as two examples of media where anisotropic wave 
turbulence is relevant. The most important results of wave turbulence are reported 
and discussed in the context of space and simulated magnetized plasmas. Important 
issues and possible spurious interpretations are eventually discussed. 



1.1 Introduction 



Wave turbulence is the stud y of the long-time stati s tical behav i or of a sea of weakly 



nonlinear dispersive waves (jZakharov et all Il992[ iNazarenkol 120111 ). The energy 



transfer between waves occurs mostly among resonant sets of waves and the resulting 
energy distribution, far from a thermodynamic equilibrium, is characterized by a 
wide power law spectrum and a high Reynolds number. This range of wavenumbers 
- the inertial range - is generally localized between large scales at which energy is 
injected in the system (sources) and small scales at which waves break or dissipate 
(sinks). Pioneering works on wave turbulence date back to the sixties when it was 
established that the stochastic initial value problem for weakly coupled wave systems 
has a natural asymptotic closure induced by the dispersiv e nature of the waves 
and the large separation of linear and nonlinear time scal es (IBennev and Saffman, 
1966 : Bennev and Newell . 1967 . 19691 ). In the meantime, Zakharov and Filonenkol 
(1966) showed that the wave kinetic equations derived from the wave turbulence 
analysis (with a Gaussian Ansatz applied to the four-point correlations of the wave 
amplitude) have exact equilibrium solutions which are the thermodynamic zero flux 
solutions but also - and more importantly - finite flux solutions which describe the 
transfer of conserved quantities b etween sources and sinks. The solution s , first 
published for isotropic turbulence ( Zakharov . 1965t Zakharov and Filonenkol 1966 ) 
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were then extended to anisotropic turbulence (jKuznetsovl 1 19721 ). 

Wave turbulence is a very common natura l phen o menon with applicati o ns, fo r 
example, in ca pillary waves 



gravity waves (IFalcon et al 



<olmakov etall booi lAbdurakhimov et all . l2008h 



su perflui d helium and pro cesses of Bose- 
Eins tein condensation iKolmakov et al. . 1995 : Lvov et al. . 2003), nonl i near op- 
tics JPvachenko et all Il992h. inertial"waves |Galtieii l2003[ iMorize et all 120051 ) or 
Alfven waves ([Galtier et all 2000; K uznetsovl |2001[ IChandranl 12005) . The most im- 



portant difference between plasmas and incompressible neutral fluids is the plethora 
of linear waves supported by the former. The direct consequence is that in weakly 
nonlinear plasmas the fundamental entities are wave s rather than the eddies of 
strong turbulence ( Kolmogorov , 1941 ; Krommesl 2002 ). In the situation discussed 
here - magnetized plasmas seen as fluids - wave and strong turbulence may coex- 
ist and therefore both waves and eddies have in fact an impact on the nonlinear 
dynamics which is strongly anisotropic. Anisotropic turbulence is particularly well 
observed in space plasmas since a magnetic field is often present on the largest scale 
of the system, like in the inner interplanetary medium w here the magnetic field 
lines form an Archimedean spiral near the equatorial plane ( Goldstein and Roberts! . 
19991 ) . at the solar surface w here coronal loops and open magnetic flux tubes are 
found IjCranmer et al.[ 12007) or in planeta r y mag netospheres where shocks and dis- 



continuities are measured ( Sahraoui et all 1200 



In the present chapter, a review is made on wave turbulence in magnetized plas- 
mas. A plasma is a gaseous state of matter in which the atoms or molecules are 
strongly ionized. Mutual electromagnetic forces between the ions and the free elec- 
trons are then playing dominant roles which adds to the complexity as compared to 
the situation in neutral atomic or molecular gases. In the situation discussed here 
the plasmas are described by the magnetohydrodynamics (MHD) approximation in 
the incompressible or compressible case. The role played by the Hall term is dis- 
cussed through the Hall MHD description as well as its small-scale limit of electron 
MHD. 

The structure of the chapter is as follows. Physical motivations for developing 
wave turbulence theories are given in Section [1.21 where we first describe multiscale 
solar wind turbulence, and then present the coronal heating problem. Section 11.31 
emphasizes the differences between strong and wave turbulence in MHD and the 
path followed historically by researches to finally obtain the MHD wave turbulence 
theories. In Section 11.41 the wave turbulence formalism is exposed with the basic 
ideas to derive the wave kinetic equations. Section 11.51 deals with the results ob- 
tained under different approximations (MHD, Hall MHD and electron MHD) in the 
incompressible or compressible case. Finally we conclude with a general discussion 
in the last Section. 
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1.2 Waves and turbulence in space plasmas 



Waves and turbulence are ubiquitous in astrophysical plasmas . Their signatures 
are found in the E arth's magnetos phere ( Sahraoui et all 200a), the solar corona 
( Chae et aL , 1998), the solar wind (Bruno and Carbone , 20051) o r the interstellar 
medium ( Elmegreen and Scald 2004 ; Scalo and Elmegreen , 2004 ). These regions 
are characterized by extremely large (magnetic) Reynolds numbers, up to 10 13 , 
with a range of available scales from 10 18 m to a few meters. 



1.2.1 Interplanetary medium 

Extensive investigations are made in the interplanetary medium (and in the Earth's 
magnetosphere which is not the subject discussed here) where many in situ space- 
crafts measurements are available. The solar wind plasma is found to be in a highly 
turbulent state with m agnetic and velocity fluctuations detected from 10~ 6 Hz up 



Smith et a 



to several hundred Hz ( Colemanl . |1968; Ro berts et all [1987; Lea mon et al 



199a 



2006). The turbulent state of the solar wind was first suggested in 1968 



( Colemanl . [1968J) when a power law behavior was reported for energy spectra with 
spectral indices lying between —1 and —2 (with the use of the Taylor "frozen- in flow" 
hypothesis). More precise measurements revealed that the spectral index at low fre- 
quency (< 1Hz) is often about —1.7 which is closer to the Kolmogorov prediction 
(IKolmogorovlll941l) fo r neut ral fluids (—5/3) rather than the Iroshnikov-Kraichnan 



prediction (|lroshnikovlll964tlKraichnanlll965) for magnetized fluids (—3/2). Alfven 
waves are also well observed since 1971 (jBelcher and Davisl . 119711 ) with a strong 
domination of antisunward propagative waves at short heliocentric distances (less 
than 1 AU). Since pure (plane) Al fven waves are e xact solutions of the ideal incom- 
pressible MHD equations (see e.g.. |Pouauetl . ll993l) . nonlinear interactions should be 
suppressed if only one type of waves is present. Therefore sunward Alfven waves, 
although subdominant, play an important role in the internal solar wind dynamics. 

The variance analysis of the magnetic field components and of its magnitude 
shows clearly that the magnetic field vector of the (polar) s olar wind has a vary - 



ing direction but with only a weak variation in magnitude (jForsvth et all 1199 



6). 



Typical values give a normalized variance of the field magnitude smaller than 10% 
whereas for the components it can be as large as 50%. In these respects, the inner 
interplanetary magnetic field may be seen as a vector lying approxim ately around a n 
Archimedean spiral direction with only weak magnitude variations ( Barnesl 1981 ). 
Solar wind anisotropy with more power perpendicular to the mean m a gnetic field 



than that parallel to it, is pointed out by data analysis ([Klein et all 119931 ) that 



provides a ratio of power up to 30. From single-point spacecraft measurements it 
is not possible to specify the exact three-dimensional form of the spectral tensor 
of the magnetic or velocity fluctuations. However, it is possible to show that the 
spacecraft-frame spectrum may depend on the angle between the local magnetic 
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field direction and the flow direction (Ho rburv et al 



20081). In the absen ce of such 



data, a quasi two-dimensional model was proposed ( Bieber et al. , 19961) in which 
wave vectors are nearly perpendicular to the large-scale magnetic field. It is found 
that about 85% of solar wind turbulence possesses a dominant 2D component. Addi- 
tionally, solar wind anisotropies is detected through radio wave scintillations which 
reveal that density spectra close to the Sun are hig hly anisotropic w i th irr egulari- 
ties stretched out mainly along the radial direction ( Armstrong et al. . 1990l) . More 
recently 

For frequencies larger than 1Hz, a steepening of the magnetic fluctuation 



power law spectra is observed over more than two decades ( Coroniti et al. . 


1982; 


Denskat et al., 


1983; 


Leamon et al. 


1998; 


Bale et al., 


2005; 


Smith et al., 2006 


i with 



a spectral index close to —2.5. This new inertial range seems to be characterized 
by a bias of the polarization suggesting that th ese fluctuations are lik ely to be 

Various 



1994) 



right-hand polarized, outward propagating waves ((Goldstein et al 
indirect lines of evidence indicate that these waves propagate at large angles to the 
background magnetic field and that the power in fluctuations parallel to the back 



ground magnetic fiel d is still less than the perpendicula r one (ICoroniti et allll982 ; 



Leamon et al. , 19981) . For these reasons, it is thought ( Stawicki et al. , 2001 ) that 



Alfven - left circularly polarized - fluctuations are suppressed by proton cyclotron 
damping and that the high frequency power law spectra are likely to consist of 
whistler waves. This scenario is supported by multi-dimensional direct numerical 
simulations of comp ressible Hall MHD turbulence in the presence of an ambient field 
rtGhosh et allll996l) where a steepening of the spectra was found on a narrow range 
of wavenumbers, and associated with the appearance of right circularly polarized 
fluctuations. This result has been confirmed numerically with a turbulent cascade 
(shell) model based on 3D Hall MHD in which a well extend ed steeper power law 
spect rum was found at scale smaller than the ion skin depth (jGaltier and Buchlin . 
20071 ). (Note that in this cascade model no mean magnetic field is assumed.) How- 
ever, the exact origin o f the change of statistical behavior is still under debate 
( Markovskii et al. . 20081 ) : for e xample, an origin from co mpressible effects is possi- 
ble in the context of Hall MHD (jAlexandrova et all 120081 ) ; a gyrokinet ic description 



is als o proposed in which kinetic Alfven waves play a central role (jHowes et al 
2008h . 



The solar wind plasma is currently the subject of a new extensive research around 
the origin of the spectral break obser yed in the magnetic fluctuations accompanied 



by an the absence of intermittency ( Kivani et al. , 20091) . We will see that wave 



turbulence may have a central role in the sense that it is a useful point of departure 
for understanding the detailed physics of solar wind turbulence. In particular, it 
gives strong results in regards to the possible multiscale behavior of magnetized 
plasmas as well as the intensity of the anisotropic transfer between modes. 
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1.2.2 Solar atmosphere 

Although it is not easy to measure directly the coronal magnetic field, it is now 
commonly acc epted that the structure o f the low solar corona is mainly due to the 
magnetic field (jAschwanden et all 120011) . The high level of coronal activity appears 
through a perpetual impulsive reorganization of magnetic structures over a large 
range of scales, from about 10 5 km until the limit of resolution, about one arcsec 
(< 726 km). The origin of the coronal reorganization is currently widely studied in 
solar physics. Information about the solar corona comes from spacecraft missions 
like SoHO/ESA or TRACE/NASA launched in the 1990s, or from more recent 
spacecrafts STEREO /NASA, Hinode/JAXA or SDO/NASA (see Fig. d). The 
most recent observations reveal that coronal loops are not yet resolved transversely 
and have to be seen as tubes made of a set of strands which radiate alternatively. 
In fact, it is very likely th at structures a t much smaller scales exist but have not 
yet been detected (see e.g., Warren . 20061) . 

Observations in UV and X-ray show a solar corona extremely hot with tempera- 
tures exceeding 10 6 K - close to hundred times the solar surface temperature. These 
coronal temperatures are highly inhomogeneous: in the quiet corona much of the 
plasma lies near 1-2 x 10 6 K and 1-8 x 10 6 K in active regions. Then, one of the 
major questions in solar physics concerns the origin of such high values of coronal 
temperature. The energy available in the photosphere - characteri zed by granules 
- is clearly sufficient to supply the total coronal losses (jPriestl |l982) which is about 



10 4 Jm~ s _1 for active regions and about one or two orders of magnitude smaller 
for the quiet corona and coronal holes where open magnetic field lines emerge. The 
main issue is thus to understand how the available photospheric energy is transferred 
and accumulated in the solar corona, and by what processes it is dissipated. 

In active region loops, analyses ma de by spectrometers show that the plasma 
velocity can reach values up to 50 km/s (jBrekke et all 11997) . The highly dynamical 
nature of some coronal loops is also pointed out by non-thermal velo cities reach- 



ing s ometimes 50 km/s as it was revealed for example by SoHO/ESA (|Chae et al 
19981 ). These observations give also evidences that the line broadening is due to 



motions which are still not resolved neither in space, with scales smaller than the 
diameter of coronal loops, nor in time, with timescales shorter than the exposure 
time of the order of few seconds. These velocity measurements are very often inter- 
preted as a signature of MHD turbulence where small scales are produced naturally 
via a nonlinear cascade of energy. In the light of the most recent observations, 
it seems fundamental to study, both theoretically and numerically, the impact of 
small-scale phenomena on the coronal heating. Note that recent Hinode/JAXA and 
SDO/NASA pictures seem to show a magnetic field controlled by plasma turbulence 
at all scales in which Alfven waves are omnipresent (see e.g., Doschek et all 2007 ; 
Nishizuka et all l2008t ICargill and De Moortel l201ll) . Thus, the turbulent activity 
of the corona is one of the key issues to understand the heating processes. 

In the framework of turbulence, the energy supplied by the photospheric motions 
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SDO/AIA 171 201 1-03-07 19:59:50 UT 



Fig. 1.1 Solar corona seen at 17.1nm by the AIA instrument onboard SDO (Solar Dynamics 
Observatory) - NASA credits. The spatial resolution of the picture is about 700km. 



and transported by Alfven waves through the corona is transferred towards smaller 
and smaller scales by nonlinear coupling between modes (the so-called energy cas- 
cade) until dissipative scales are reached from which the energy is converted into 
heating. The main coronal structures considered in such a scenario are the magnetic 
loops which cover the solar surface. Each loop is basically an anisotropic bipolar 
structure anchored in the photosphere. It forms a tube of magnetic fields in which 
the dense and hot matter is confined. Because a strong guiding magnetic field (Bo) 
is present, the nonlinear cascade that occurs is strongly anisotropic with small scales 
mainly developed in the Bo transverse pl anes. Most of the models publis hed deals 
with isotropic MHD turbulence (see e.g.. iHendrix and Van HovenL Il996l ) and it is 
only very recently that a nisotropy has been included in turbulent heating models 
(jBuchlin and Vellil . l200i . 
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The latest observations show that waves and turbulence are among the main 
ingredients of the solar coronal activity. Weak MHD turbulence is now invoked has 
a possible regime for some coronal loops since a very s mall ratio is found between th e 
fluctuating magnetic field and the axial component ( Rappazzo et al. . 2007 , 20081 ). 
Inspired by the observations an d by recent direct numerical simulations of 3D MHD 
~[|2008bl) . an analytical model of coronal structures has been 
2008cl ) where the heating is seen as the end product of a 



turbulence (IBigot et al 



proposed (Big ot et al 



wave turbulent cascade. Surprisingly, the heating rate found is non negligible and 
may explain the observational predictions. 

The coronal heating problem also concerns the regions where the fast solar wind 
is pr oduced, i.e. the coronal holes ([Hollweg and Isenbergj . l2002t ICranmer et al. . 
20071 ). Observations seem to show that the heating affects preferentially the ions 



in the direction perpendicular to the mean magnetic field. The electrons are much 
cool er than the ions, w ith temperatures generally less than or close to 10 6 K (see 
e.q.. lDavid et all Il99a ). Additionally, the heavy ions become hotter than the pro- 



tons within a solar radius of the coronal base. Ion cyclotron waves could be the 
agent which heats the coronal ions and accelerates the fast wind. Naturally the 
question of the origin of these high frequency waves arises. Among different sce- 
narios, turbulence appears to be a natural and efficient mechanism to produce ion 
cyclotron waves. In this case, the Alfven waves launched at low altitude with fre- 
quencies in the MHD range, would develop a turbulent cascade to finally degenerate 
and produce ion cyclotron waves at much higher frequencies. In that context, the 
wave turbulence regime was considered in the weakly compressible MHD case at 
low-/3 plasmas (where /3 is the ratio between the thermal and magnetic pressure) 
in orde r to analyze the n onlinear three-wave interaction transfer to high frequency 
waves dChandranl . 120051) . The wave turbulence calculation shows - in absence of 
slow magnetosonic waves - that MHD turbulence provides a convincing explanation 
for the anisotropic ion heating. 



1.3 Turbulence and anisotropy 

This Section is first devoted to the comparison between wave and strong turbulence. 
In particular, we will see how the theoretical questions addressed at the end of the 
20th century have led to the emergence of a large number of papers on wave turbu- 
lence in magnetized plasmas and to many efforts to characterize the fundamental 
role of anisotropy. 



1.3.1 Navier-Stokes turbulence 

Navier-Stokes turbulence is basically a strong turbulence problem in which it is 
impossible to perform a (non trivial and) consistent linearization of the equa- 
tions against a stationary homogeneous background. We remind that wave tur- 
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bulence demands the existence of linear (dispersive) propagative waves as well as 
a large separation of line ar and nonlinear (eddy-turnover) ti me scales (see e.g. , 
Bennev and Newell . 1969 ). In his third 1941 turbulence paper, Kolmogorov ( 1941 ) 
found that an exact and nontrivial relation may be derived from Navier - Stoke s 
equations for the third-order longitudinal structure function (|Kolmogorovl Il94lh : 
it is the so-called 4/5s law. Because of the rarity of such results, the Kolmogoroy's 
4 /5s law is considered as one of the most important results in turbulence ([Frisch . 
19951 ). Basically, this law makes the following link in the 3D physical space between 
a two-point measurement, separated by a distance r, and the distance itself] 



-e r ■■ 



Vh) 



(1.1) 



where () denotes an ensemble average, the longitudinal direction L is the one along 
the vector separation r, v is the velocity and e v is the mean (kinetic) energy injection 
or dissipation rate per unit mass. To o btain this exact relation, the assumptions 
of homogeneity and isotropy are made ( Batchelor . 19531 ). The former assumption 
is satisfied as long as we are at the heart of the fluid (far from the boundaries) 
and the latter is also satisfied if no external agent (like, for example, rotation or 
stratification) is present. Additionally, the long time limit is considered for which a 
stationary state is reached with a finite e v and the infinite Reynolds number limit 
(i.e. the viscosity v — >■ 0) is taken for which the mean energy dissipation rate per 
unit mass tends to a finite positive limit. Therefore, the exact prediction is valid in 
the asymptotic limit of a large inert ial range. Th e Kolmogorov law is well supported 
by the experimental data (see e. 



FrischL 11995 ) 



The 4/5s law is a fundamental result used to develop heuristic spectral scaling 
laws like the famous - but not exact - 5/3-Kolmogorov energy spectrum. This 
point makes a fundamental difference with wave turbulence where the power law 
spectra found are exact solutions of the asymptotically exact wave turbulence equa- 
tions. Nevertheless, the term "Kolmogorov theory" is often associated to the —5/3 
spectrum since there exists a theory behind in the physical space. 



1.3.2 Incompressible MHD Turbulence 
1.3.2.1 Strong turbulence 

The wave turbulence regime exists in incompressible MHD. The main reason is 
that Alfven waves are linear solutions when a stationary homogeneous background 
magnetic field Bo is applied. This statement seems to be obvious but we will see 
that the problem is subtle and the existence of an Alfven wave turbulence theory 
was the subject of many discussions basically because those waves are only pseudo- 
dispersive (i.e. , the frequency lo is proportional to a wavenumber). 



lr The Kolmogorov 4/5 s law can be written in an equivalent way as a 4/3s law l lAntonia et all 
Il997l : iRasmussenl. Il999t) . namely: -(4/3)e"r = {(v' L - v^'EiH ~ v i) 2 )- 
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The question of the existence of an exact relation for third order structure func- 
tions is naturally addressed for strong (without Bo) MHD t urbulence. The answer 
was given by Politano & Pouquet only in 1998 (see also, IChandrasekhan, 119511 ) 
for incompressible MHD turbulence. The presence of the magnetic field and its 
coupling with the velocity field renders the problem more difficult and, in practice, 
we are dealing with a couple of equations. In this case, the possible formulation in 
3D is the 4/3's law 

4 



-e±r = {{z'T - zDYuf - Z ff) , (1.2) 



where the direction L is still the one along the vector separation r, = v ± b 
is the Elsasser fields (with b normalized to a velocity field) and e is the mean 
energy dissipation rate per unit mass associated to the Elsasser energies. To obtain 
these exact relations, the assumptions of homogeneity and isotropy are still made, 
and the long time limit for which a stationary state is reached with a finite e is 
also considered. The infinite kinetic and magnetic Reynolds number limit (y — > 
and the magnetic diffusivity r\ —¥ 0) for which the mean energy dissipation rates 
per unit mass have a finite positive limit is eventually taken. Therefore, the exact 
prediction is again valid, at first order, in a wide inertial range. This prediction is 
currently often used in the literature to analyze space plasma data and det e rmine 
the local heating rate e: ± in the solar wind (see e.g., ISorriso-Valvo et all 12007 : 
MacBride et~ail 120081 ) . 



1.3.2.2 Iroshnikov-Kraichnan spectrum 

The isotropy assumption used to derive the 4/3's law is stronger for magnetized 
than neutral fluids since in most of the situations encountered in astrophysics a 
large scale magnetic field Bo is present which leads to anisotropy (see Section [172} . 
Although this law is a fundamental result that may be used to develop a heuristic 
spectral scaling law, the role of Bq has to be clarified. Indeed, we have now two 
time-scales: the eddy-turnover time and the Alfven time. The former is similar to 
the eddy-turnover time in Navier-Stokes turbulence and may be associated to the 
distortion of wave packets (the basic entity in MHD), whereas the latter may be 
seen as the duration of interaction between two counter-propagating Alfven wave 
packets (see Fig |1.2l) . During a collision, there is a deformation of the wave packets 
in such a way that energy is transferred mainly at smaller scales. The multiplicity 
of collisions leads to the formation of a well extended power law energy spectrum 
whose index lies between —5/3 (Kolmogorov prediction) and —3/2 (Iroshnikov- 
Kraichnan prediction) according to the phenomenology used, i.e. with or without 
the Alfven wave effect. Note the slightly different approaches followed by Iroshnikov 
and Kraichnan to derive the —3/2 spectrum. In the former case the presence of 
a strong magnetic field is explicitly assumed whereas it is not in the latter case 
where it is claimed that the small-scale fluctuations see the large-scales - in the 
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sub-inertial range - as a spatially uniform magnetic field. However, they both 
assumed isotropy to finally proposed the so-called Iroshnikov-Kraichnan spectrum 
for MHD turbulence. It is important to note that the exact isotropic relation (|1.2j) 
in physical space corresponds dimensionally to a —5/3 energy spectrum since it is 
cubic in z (or in v and b) and linear in r. It is therefore less justified to use the term 
"theory" for the Iroshnikov-Kraichnan spectrum than for the Kolmogorov one for 
neutral fluids. 



1.3.2.3 Breakdown of isotropy 

The weakness of the Iroshnikov-Kraichnan phenomenology is the apparent contra- 
diction between the (in) direct presence of a strong uniform magnetic field and the 
assumption of isotropy. An important difference between neutral and magnetized 
fluids is the impossibility in the latter case to remove a large-scale (magnetic) field 
by a Galilean transform. The role of a uniform magnetic field has been widely 



~- ~ -~ ~- , 

e.o.. Montgomery and Turner. 1981: 


— ~ — , — ~~ 

Shcbalin et al.. 


Ne and BhattacharieeL 1996: 


Verma . 


2004). At stron 



clearly established results is the bi-dimensionalization of MHD turbulent flows with 
a strong reduction of nonlinear transfers along Bo- The consequence is an energy 
concentration near the plane k ■ Bq = 0, a result illust rated later on by dir ect 



1983) 



numerical simulations in two and three space dimensions (phebalin et al. 

The effects of a strong uniform magnetic field may be handled through an anal- 
ysis of resonant triadic interactions (|Shebalin et all 119831 ) between the wavevectors 
(k, p, q) which satisfy the vectorial relation 



k = p + q, 

whereas the associated wave frequencies satisfy the scalar relation 

u>(k) = w(p) + w(q) . 
For incompressible MHD, the Alfven frequency is 

w(k) = ±k-B =±k\\B , 



(1.3) 
(1.4) 
(1.5) 




Fig. 1.2 Alfven wave packets propagating along a large scale magnetic field Bq. 
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where || defines the direction along Bo (_L will be the perpendicular direction to Bo 
which is written in velocity unit). The solution of the three- wave resonance con- 
dition gives for example, qu = 0, which implies a spectral transfer only in the per- 
pendicular direction. For a strength of Bq well above the r.m.s. level of the kinetic 
and magnetic fluctuations, the nonlinear interactions of Alfven wave packets may 
dominate the dynamics of the MHD flow leading to the regime of wave turbulence 
(see Section 11.41) where the energy transfer, stemming from three-wave resonant 
interactions, can only increase the perpendicular component of the wavevectors, 
while the nonlinear transfers are completely inhibited along Bo- The end result is 
a strongly anisotropic flow. 



1.3.2.4 Emergence of anisotropic laws 

An important issue discussed in the literature is the relationship be tween perpen- 
dicular and parallel scales in anisotrop i c MHD turbulence ( see e.g.. iHigdonl . Il984 : 



Goldreich and Sridhai 



1995 



Bol dvrev 



2006 



Sridhar, 2010). In order to take into 



account anisotropy, iGoldreich and Sridharl (|1995l ) proposed a heuristic model based 
on the conjecture of a critical balance between the Alfven and the eddy-turnover 
times, which are respectively 



t a ~£\\/B 



and 



Teddy 



(1.6) 



(1.7) 



where H\\ and £j_ are typical length scales parallel and perpendicular to Bo. The 
conjecture says that in the inertial range we have, ta = T e ddy The latter relation 
leads trivially to, ue ~ BqI_l/£\\. Following the Kolmogorov arguments, one ends 
up trivially with the Kolmogorov energy spectrum 



E(kj_,k 



-5/3 



(where k = (kj_, fc||), k±_ = |kj_|, / J E{k^,k\\)dk^dk\\ — J E(kj_)dk_i 
the non trivial anisotropic scaling law (with uJ/T e ddy = constant) 



(1.8) 



and with 



II 



,2/3 



(1.9) 



This heuristic prediction means that anisotropy is stronger at smaller s cales. 

A generalization of this result was proposed by iGaltier et al.1 (|2005T ) in order to 
model MHD flows both in the wave and strong turbulence regimes, as well as for 
the transition between them. In this heuristic model, the constant time-scale ratio 
X = ta/ T e ddy is not assumed to be necessarily equal to unity. The relaxation of the 
constraint (x = 1) allows us to recover the anisotropic scaling law (jl.9p which now 
includes Bq, 

k^k^/Bo, (1.10) 
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and to find a universal prediction for the total energy spectrum 



E(k±,k\\) ~ fcj"fc 



(1.11) 



with 



(1.12) 

^~eddy /X 



3a + 2(3 = 7 . 

Note that a classical calculation with a transfer time r = T^ ddy /TA = 
leads trivially to the spectrum (11.81) : the choice of r fixes irreversibly the spec- 
trum. It is only when the ansatz (jl.lll) is introduced that the non trivial re- 
lation (I1.12p emerges; this ansatz is believed to be weak since power laws are 



see 



e.g. 



Cho and Vishniac 



Shaikh and Zank . 



2007 



2000; 



Bigot et all . l2008bF . the anisotiopic scaling law between 





2001: 


Ne et al. 





2003; 



parallel and perpendicular scales (I1.9P seems to be a robust result and an approx- 
imately constant ratio Xj generally smaller than one, is found between the Alfven 
and the eddy-turnover times. This sub-critical value of x implies therefore a dy- 
namics mainly driven by Alfven waves interactions. Note that the presence of Bq 
in relation (|1.10j) shows the convergence towards wave turbulence (Bq — > +oo, with 
respect to the fluctuations) for which the parallel transfer is totally frozen. 

The ques tion of the spectral indices i s still a challenging problem in anisotropic 
turbule nce (jSagaut and Camborl 2008 ). The main conclusion of Bigot et al 



(|2008bl ) is that the difficulty to make the measurements is generally underestimated 
in a sense that the scaling prediction in k± may change significantly when E(k±, ) 
is plotted at a given k\\ instead of E(k±): indeed, the slow mode E(k±, k\\ = 0) may 
play a singular role in the dynamics with a scaling in k±_ different from the one given 
by the 3D modes E(k_\_,k\\ > 0). This comment holds primarily for direct numer- 
ical simulations where technically it is currently possible to make this distinction, 
observations being still far from this possibility. This point will be further discussed 
in the last Section. Note finally that all these spectral predictions suffer from rig- 
orous justifications and the word "theory" that we find very often in the literature 
is not justified at all. A breakthrough could be achieved if one could develop the 
equivalent of an exact 4/5's law for anisotropic MHD turbulence. Then, a dimen- 
sional derivation from this law could lead to an anisotropic spectral prediction and 
in some sense, for the first time, the possibility of having a theoretical link between 
strong and wave turbulence. Recently, an attempt has been made in that d irection 
by using the idea of critical balance for third-order moments (|Galtierl . 120121 ) . More- 
over, it was shown that the introduction of the dynamic alignment conje cture into 
the exact relation for third-order momen ts ( Politano and Pouauet . 199Sl) may give 
a kj 3 ^ 2 spectrum ( Boldvrev et al . 2009f) . 



I. 3.3 Towards an Alfven wave turbulence theory 

In vi ew of the importance o f aniso tropy in natural magnetized plasma (see Section 

II. 2p , ISridhar and Goldreichl (jl994l ) suggested that a plasma evolving in a medium 
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permeated by a strong uniform magnetic field and in the regime of wave turbulence 
is characterized by four-wave nonlinear interactions. The essence of wave turbu- 
lence is the statistical study of large ensembles of weakly interacting waves via a 
systematic asymptotic expansion in powers of small nonlinearity. This technique 
leads finally to the exact derivation of wave kinetic equations for the invariants of 
the system like the energy spectrum (see Section \T^i . In MHD, it is the presence 
of a strong uniform magnetic field Bo that allows to introduce a small parame- 
ter in the system, namely the ratio between the fluctuating fields and Bq. The 
result found by Sridhar and Goldreich in 1994 implies that the asymptotic devel- 
opment has no solution at the main order and that we need to go to the next 
order to describe Alfven wave turbulence. Several a r ticles , using a phenomenology 
(see e.g., Montgomery and Matthaeusl . 1995 ; Verma , 2004 ) or a rigorous treatment 
(|Ng and Bhattacharieei Il996l ). were published to contest this conclusion and sus- 
tain the non trivial character of the three-wave interactions . In response, a de- 



tailed theory was fi nally given in 2000 (jGaltier et all 120001 ; iNazarenko et all 12001 



Galtier et all 120021) whose main prediction may be derived heuristically in few lines 
as follows. According to Fig. 11.21 the main process which happens in Alfven wave 
turbulence is the stochastic collisions of wavepackets. To find the transfer time and 
then the energy spectrum, first we shall evaluate the modification of a wavepacket 
produced by one collision. We have (for simplicity we only consider the balanced 



case for which z 



± 



- z e — u e ~ bi) 
zt{t + t a ) ~ z e (t) 



T A —~ Ze (t)+T A -. 



(1.13) 



where t a is the duration of one collision; in other words, after one collision the 
distortion of a wavepacket is AiZ£ <~ t a z^/£±. This distortion is going to increase 
with time in such a way that after N stochastic collisions the cumulative effect may 
be evaluated like a random walk 



N 

i=i 



T A - 



(1.14) 



The transfer time Tt r that we are looking for is the one for which the cumulative 
distortion is of the order of one, i.e. of the order of the wavepacket itself 



Zt ~ T A 



z l / T tr 



e± V r A ' 



(1.15) 



then we obtain 



T A z] 



eddy 
T A 



(1.16) 



A classical calculation, with e ~ z|/r tI ., leads finally to the energy spectrum 



(1.17) 
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where ku has to be seen as a parameter since no transfer along the parallel direction 
is expected (see Section fl .3.2.3^ . Note that this demonstration is the one tradi- 
tionally used for deriving the Iroshnikov-Kraichnan spectrum, but since in this case 
isotropy is assumed a E(k) ~ ^eBq k~ 3 ^ 2 is predicted. 

First signatures that may be attributed to wave turbulence were found by 
Perez and Boldvrevl ( 2008^) in numerical simulations of a reduced form of the MHD 
equations (jGaltier and Chandran , 2006 ). The detection of the wave turbulence 
regime from direct numerical simulations of MHD equations is still a difficult task 
but recent results have been obtained in wh ic h temporal, structural an d spectral 



signatures are reported ( Bigot et al. , 2008al lbl: Bigot and Galtier , 2011 ). Current 



efforts are also made to analyze the effects of oth er inviscid invariants, like the 
cross-helic i ty, on the scaling laws of wave turbulence ( Lithwick and Goldreich , 20031 : 
Chandranl 120081 ). It is worth noting that these works on imbalanced wave turbu- 
lence have been followed by (sometimes controversial) research investigations in the 
strong turbulence regime which is of great relevance for the solar wind (see e.g . 
Chandran et al. , 2009 : Beresnvak and Lazarianl 2010l : Podesta and Bhattachariee . 
2010). 



1.3.4 Wave turbulence in compressible MHD 

Most of the investigations devoted to wave turbulence refers to isotropic media 
where the well known conformal transform proposed by Zakharov and Filonenkol 
(1966) may be applied to find the so-called Kolmogorov- Zakharov spectra (see 
Section 11.4.41) . (Surprisingly a similar transform was used in the meantime by 



Kraichnan ( 1967t) to investigate the problem of 2D turbulence.) The introduction 
of anisotropy in plasmas was stu died to a smaller e xtent: the first example is given 
by magnetized ion-so und waves (IKuznetsovl . I1972T) . The compressible MHD case 
was analyzed later by IKuznetsovl (|200ll ) for a situation where the plasma (thermal) 
pressure is small compared with the magnetic pressure (small ft limit). In this 
case, the main nonlinear interaction involving MHD waves is the scattering of a 
fast magneto- acoustic and Alfven waves on slow magneto-acoustic waves. In this 
process, the fast and Alfven waves act as high-frequency waves with respect to the 
slow waves. (To simplify the analysis other three-wave interaction processes that 
do not involve slow waves are neglected.) 

A variant of the wave turbulence analysis in compressible MHD was proposed 
bv lChandranl l|2005[) in the limit of a small /3 in which the slow waves are neglected 
and a constant density is imposed. The other (mathematical) difference is that 
the Hamiltonian formalism was used in the former analysis whereas an Eulerian 
description was employed in the latter case. Because the compressible regime is 
much more difficult to analyze than the incompressible one, simplifications have 
been made and to date, no general theory has been proposed for compressible MHD 
wave turbulence. 
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1.3.5 Wave turbulence in Hall and electron MHD 

Modeling the physics of a plasma beyond the MHD approximation, namely on spa- 
tial scales shorter than the ion inertial length di (but larger than the electron inertial 
length d e ) and time scales shorter than the ion cyclotron period l/w c j, is a highly 
challenging problem even in the fluid approximation. (For kinetic models we can 
mention the gyrokinetic approximation useful for weakly collisio nal plasmas in which 
time scales are supposed to be much larger than l/ui C i (see e.g., Schckochihi n et al 



20091 ).) In that context the electron MHD approximation ([Kingsep et all I1990T ) is 



often used: in such a limit one assumes that ions do not have time to follow elec- 
trons and provide a static homogeneous background on which electrons move. The 
electron MHD approximation is particularly relevant in the context of collision- 
less magnetic reconnection where a diffusion region is developed with multiscale 
structures corresponding to ion and electron characteristic lengths (jBiskampl . 11997 ; 
Bhattacharieei |2004 lYamadal . l2007h . 

An important issue i n electron MHD turbu lence is about the impact of whistler 
waves on the dynamics. iBiskamp et al.l (|199.9) argued that although whistler wave 
propagation effects are non negligible in electron MHD turbulence, the local spectral 
energy transfer process is independent of the linear wave dispersion relation and the 
energy spectrum may be predicted by a Kolmogorov type argument. Direct numer- 
ical simulati ons were used to illustrate the theory b ut no mean magnetic field was 



introduced ( Biskamp et al. . 19961 : Ng et al. . 2003 ). Dastgeer et al. ( 2000 ) investi- 



gated the turbulence regime in the presence of a moderate background magnetic 
field Bq and provided convincing numerical evidence that turbulence is anisotropic. 
It was argued that although whistler waves may appear to play a negligible role 
in determining the spectral index, they are important in setting up an anisotropic 
turbulent c ascade. The whistler wave turbulence re gime was then investigated theo- 
retically by Galtier and Bhattachariee ( 20031 2005 ) in the limit Bq +oo (with re- 
spect to the fluctuations). It was shown that similarly to the MHD case, anisotropy 
is a central feature of such a turbulence. Att e mpts to find a n anisotropic law for 
electron MHD was made bv lCho and Vishniad (|2004l ) (see also. lGaltier et all 120051 ) 

1/3 

and a scaling relation in fc|| ~ was found both heuristically and numerically 

Hall MHD is an extension of the standard MHD where the ion inertia is retained 
in Ohm's law. It provides a multiscale description of magnetized plasmas from which 
both standard and electron MHD approximations may be recovered. Hall MHD is 
often u sed to understand, for example, th e magnetic field evolution in neutron star 



crust s (jGoldreich and Reiseneggerl Il992l). the turbulent dynamo (jMininni et al 



20031 ). the formation of f ilaments (Laveder. Passot and Sulcm, 2002) , the multiscale 
solar wind turbulence ( Ghosh et al. . 19961: Krishan and Mahaian . 2004 : Galtier , 
2006allbh . or the dynamics of the magnetosheath ([Belmont and Rezeaul 120011) . 
Anisotropy in Hall MHD is clearly less understood than in MHD or electron MHD 
mainly because the numerical treatment is more limited since a wide range of scales 
is necessary to detect any multiscale effects. From a theoretical point of view, 
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it is only recently that a wave turbulence theory has been achieved for the in- 
compressible case (jGaltierl l2006bl : ISahraoui et all booi . For such a turbulence, 



a global tendency towards anisotropy was found (with, however, a weaker effect 
at intermediate scales) with nonlinear transfers preferentially in the direction per- 
pendicular to the external magnetic field B . The energy spectrum is character- 
ized by two inertial ranges, separated by a knee, which are exact solutions of the 
wave kinetic equations. The position of the knee corresponds to the scale where 
the Hall term becomes sub/dominant. To date, compressible Hall MHD is still 
an open problem in the regime of w ave turbulence. (A first step was made by 
Sahraoui. Belmont and Rezeau ( 20031) who found the Hamiltonian system.) 



1.4 Wave turbulence formalism 
1.4.1 Wave amplitude equation 

Wave turbulence is the study of the long time statistical behavior of a sea of weakly 
nonlinear dispersive waves. It is described by wave kinetic equations. In this section 
we present the wave turbulence formalism which leads to these nonlinear equations. 
We shall use the inviscid model equation 

^M=£(u)+ £ AA(u,u), (1.18) 

where u is a stationary random vector, £ is a linear operator which insures that 
waves are solutions of the linear problem, and AT is a quadratic nonlinear operator 
(like for MHD- type fluids). The factor e is a small parameter (0 < e < 1) which 
will be used for the weakly nonlinear expansion. For all the applications considered 
here, the smallness of the nonlinearities is the result of the presence of a strong 
uniform magnetic field Bo; the operator £ is thus proportional to Bq. 
We introduce the 3D direct and inverse Fourier transforms 



u(x, t) = / A(k, t) exp(ik • x)c0c , 
A (M) = 77TTT / u(x,t)exp(-ik-x)dx. 



(1.19) 
(1.20) 



Therefore, a Fourier transform of equation (jl.181) gives for the j-component 

' d_ 

.dt 



|+ lW (k))^(k,t)= (1.21) 



£ / H^lA m (p,t)A n (ci,t)S(k-p-q)dpdq, 

where w(k) = ujk is given by the appropriate dispersion relation (with in general 
lo{— k) = — w(k)) and % is a symmetric function in its vector arguments which 
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basically depends on the quadratic nonlinear operator J\f. Note the use of the 
Einstein's notation. We introduce 

A(k,t) =a(k,t)e- ia;f!t , (1.22) 

and obtain in the interaction representation 

^M= s [ ^a m (p)a„(q)e in *-*r5 fc , Mt ipriq, (1.23) 
at j R6 

where the Dirac delta function §k,pq = S(k— p— q) and f2fe >pg = ^k—u p —Lu q ; the time 
dependence in fields, a, is omitted for simplicity. Relation (|1.23l) is the wave am- 
plitude equation whose dependence in e means that weak nonlinearities will modify 
only slowly in time the wave amplitude. By nature, the problems considered here 
(in MHD, electron and Hall MHD) involve mainly three- wave interaction processes 
as it is expected by the form of the wave amplitude equation. The exponentially 
oscillating term is essential for the asymptotic closure since we are interested in the 
long time statistical behavior for which the nonlinear transfer time is much greater 
than the wave period. In such a limit most of the nonlinear terms will be destroyed 
by random phase mixing and only a few of them - called the resonance terms - 
will survive. Before going to the statistical formalism, we note the following general 
properties that will be used 

«S=(«S: P "T. (1-24) 
■Hjrnn is symmetric in (p, q) and (m, n) , (1-25) 

H°™=0, (1.26) 

where, *, stands for the complex conjugate. 

1.4.2 Statistics and asymptotics 

We turn now to the statistical description, introduce the ensemble average (...) and 
define the density tensor for homogeneous turbulence 

fty(k')<5(k + k') = (oj-flcJaj/pO) • (1-27) 

We also assume that on average (u(x, t)) = which leads to the relation / Hj^ l = 0. 
From the nonlinear equation (|1.23[) . we find 

d ^ k + ^ - L^f-^m + /a,(k)^\ = (1.28) 



e 



dt \ 3 v ' dt I \ JK ' dt 

^mn< a m(P) a n(q)a,'(k'))e infc -*4, P ^pdq 



+ 



e / <^(a„ l (p)a n (q)a J (k))e in -'-*4', P9 dpdq. 
Jr. 6 
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A hierarchy of equations will clearly appear which gives for the third order moment 
equation 

a(a 3 -(k)o i /(k')o i /,(k")) 



dt 



(1.29) 



£ / H^( am (p)a„(q)a J v(k')a J »(k"))e l0fc -*4^dpdq 
{(k,j)^(k',/)}dpdq 



e^ 6 {(k",j")o(k',j')}rfpdq, 



where in the right hand side the second line means an interchange in the notations 
between two pairs with the first line as a reference, and the third line means also an 
interchange in the notations between two pairs with the second line as a reference. 
At this stage, we may write the fourth order moment in terms of a sum of the 
fourth order cumulant plus products of seco nd order ones, but a natural closure 
arises for times asymptot ically large (see e.g., Newell et al. . 2001 ; Nazarenkcl 2011 ; 



Newell and RumpJ . 120111 ). In this case, several terms do not contribute at large 
times like, in particular, the fourth order cumulant which is not a resonant term. 
In other words, the nonlinear regeneration of third order moments depends essen- 
tially on products of second order moments. The time scale separation imposes a 
con dition of applicabi lity of wave turbulence which has to be checked in fine (see 
Nazarenkol 120071 ). After integration in time, we are left with 



e.g.. 



where 



(a i (k)o,v(k')o J v,(k // )>= (1-30) 
^ / ^((a m (p)a„(q))(a / (k')a,»(k")) + (a m (p)a / (k')><a„(q)a j «(k")) 
+(a m (p)a 3 -«(k"))(a n (q)o 3 v(k'))^A(nfc j pg)5fc jP ,dpdq 
ej^ {(k,j) o (k',f)]dpdq + ej ^ {(k",j") o (k' ,j')}dpdq, 



pq I 



t»l/a 



J^k.pgt ^£ _ 



in 



k,pq 



(1.31) 



The same convention as in (|1.29p is used. After integration in wave vectors p and 
q and simplification, we get 



(a 3 -(k)o J v(k / )o J v,(k")) 



(1.32) 



eA(ilkk'k")8kk'k" 
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(u)-«-«'q mj , (k') qnj „ (k") + U)-«'-*q mj „ (k")q nj , (k') 

+< m r k ''^(k)g n/ '(k'0+<„^^ k g m ,»(k'')^(k) 

+U$-*- k 'q mj (k)« ni , (k') + H*:^-\ mj . (V)q nj (k) 
The symmetries (11.25)) lead to 

(a i (k)o,v(k')o,v,(k // )>= (1-33) 

2 £ A(fi fefe , fe »)W' (« k ™'~ k 'W(k')9ry"( k ") 

+^^ k ''^(k) 9 „ J ,,(k'o+^^ k '^•(k) f3 „ i^ (k'; 

The latter expression may be introduced into (ll.28[) . We take the long time limit 
(which introduces irreversibility) and find 

A(x) -+n5(x)+iV(l/x), (1.34) 

with V the principal value of the integral. We finally obtain the asymptotically 
exact wave kinetic equations 

(HP;r k 9r„(q)<Z /s (k) + ^- pk 9rm (p)^ s (k) + %/ r ; p - \ m ( P )? s „(q)) dpdq . 

These general 3D wave kinetic equations are valid in principle for any situation 
where three- wave interaction processes are dominant; only the form of % has to be 
adapted to the problem. Equation for the (total) energy is obtained by taking the 
trace of the tensor density, ^-(k), whereas other inviscid invariants are found by 
including non diagonal terms. 



1.4.3 Wave kinetic equations 

Equation (11.35)) is the wave kinetic equation for the spectral tensor components. We 
see that the nonlinear transfer is based on a resonance mechanism since we need to 
satisfy the relations 

LU k = Up+UJq, (1.36) 

k = p + q. (1.37) 

The solutions define the resonant manifolds which may have different forms accord- 
ing to the flow. For example in the limit of weakly compressible MHD, when the 
sound speed is much greater than the Alfven speed, it is possible to find (for the 
shape of the resonant manifolds in the 3D k-space) spheres or tilted planes for Fast- 
Fast- Alfven and Fast-Fast-Slow wave interactions, and rays (a degenerescence of the 
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resonant manifolds) for Fast-Fast-Fast wave interactions (jGaltier et all 120011 ). We 
also find planes perpendicular to the uniform magnetic field B for Slow-Slow-Slow, 
Slow-Slow- Alfven, Slow- Alfven- Alfven or Alfven- Alfven- Alfven wave interactions; it 
is a similar situation to incompressible MHD turbulence (since, at first order, slow 
waves have the same frequency as Alfven waves) for which the resonant manifolds 
foliate the Fourier space. 

The representation of the resonant manifolds is always interesting since it gives 
an idea of how the spectral densities can be redistributed along (or transverse to) the 
mean magnetic held direction whose main e ffect is the nonlinear transfer reduction 

19961) 



along its direction ([Matthaeus et al 



The previous findi ng was confirmed 



by a deta i led an alysis of wave compressible MHD turbulence (jKuznetsovl . 12001 



Chandranl 120051) in the small /3 limit where the wave kinetic equations were derived 
as well as their exact power law solutions. The situation for electron and Hall MHD 
is more subtle and there is no simple picture for the resonant manifolds like in 
MHD. In this case, it is nevertheless important to check if the resonance condition 
allows simple particular solutions in order to justify the domination of three-wave 
interaction processes over higher order (four-wave) processes. 

The form of the wave kinetic equations (jl.35|) is the most general one for a dis- 
persive problem as whistler wave turbulence in electron MHD. The incompressible 
MHD system constitutes a unique example of pseudo-dispersive waves for which 
wave turbulence applies. In this particular case, some symmetries are lost and the 
principal value terms remain present. For that reason, incompressible MHD may 
be seen as a singular limit of incompressible Hall MHD (jGaltierl . liooibl ) . 



1.4.4 Finite flux solutions 

The most spectacular result of the wave turbulence theory is its ability to pro- 
vide exact finite flux solutions. These solutions are found afte r applying to the 
wave kinetic equations a conformal transform proposed first by Zakharov ( 19651 ) 



for isotropic turbulence: it is the so-called Kolmogorov-Zakharov spectra. Because 
anisotropy is almost always present in magnetized plasm as, a bi-homogeneous con- 
formal transform is more appropriate (jKuznetsovl . Il972h . This operation can only 
be performed if first one assumes axisymmetry. With this assumption the wave 
kinetic equations (|1 .35|) write 

^ y( ^' fc||) = 4^ 2 / r4 W(«*«)«SS (i-38) 

(^m7^ k 9™ (5-L . 9H )Qj> s {k± , k± ) + H^~ pk g rm (p±,p\\)qf s {k±, k\\ ) 

J P\\ )q~sn 

where qjj'(kj_, kn) = g^y (k)/(27rfcj_) and % is a geometric operator. Note that we 
have performed an integration over the polar angle. Except for incompressible MHD 
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for which the wave kinetic equations simplify thanks to the absence of nonlinear 
transfer along the parallel direction (along Bo), in general we have to deal with a dy- 
namics in the perpendicular and parallel directions with a relatively higher transfer 
transverse to Bo than along it. In this case, we may write the wave kinetic equa- 
tions in the limit k± ^> kn . The system of integro-differential equations obtained is 
then sufficiently reduced to allow us to extract the exact power law solutions. We 
perform the conformal transformation to the equations for the invariant spectral 
densities (total energy, magnetic helicity...) 



P± -> k ±/P±- > 
qj_ -t k ± q ± /p ± 

P\\ -> fejj/i?|| . 

q\\ -> k \\i\\/p\\ ' 



(1.39) 



and we search for stationary solutions in the power law form fcj n rvjj~ m . Basically two 
types of solutions are found: the fluxless solution, also called the thermodynamic 
equilibrium solution, which corresponds to the equipartition state for which the 
flux is zero, and the finite flux solution which is the most interesting one. During 
the last decades many papers have been devoted to the find ing of these finite flux 



solutions for isotr opic as well as anisotropic wave turbulence (jZakharov et al 
Nazarenkolboilh . 



1991 



Recently, and thanks to high numerical resolutions, a new challenge has ap- 
peared in wave turb ulence: the study of incompressible Alfven wave turbulence 
(jGaltier et all 120001 ) reveals that the development of the finite flux solution of the 
wave kinetic equations is preceded by a front characterized by a significantly steeper 
scaling law. An illustration is given in Fig. ll.3l for the balance case [E = E + = E~): 
the temporal evolution of the energy spectrum reveals the formation of a k ± front 
before the establishment of the Kolmogorov-Z akharov solution in kj 2 . This find ing 
is in contradiction with the theory proposed bv lFalkovich and Shafarenkol (|199ll) on 
the nonlinear front propagation where the authors claimed that the Kolmogorov- 
Zakharov spectrum should be formed right behind the propagating front. The same 
obser vation was also mad e for the inverse cascade in the nonline ar Schrodinger equa- 



tion ( Lacaze et all 120011 ) and a detailed analysis was given by IConnaughton et al 



(2003) in the simplified case of strongly local interactio ns for which the usual way e 
kinetic equations become simple PDEs (for MHD, see. lGaltier and Buchlin . 2010). 
In spite of these attempts, the anomalous scaling is still an open problem without 
clear physical and proper mathematical answer. 
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k, 



Fig. 1.3 Time propagation (from left to right) of the compensated energy spectrum for incom- 
pressible Alfven wave turbulence (E = B+ = E~): an anomalous scaling is observed before the 
formation of the Kolmogorov-Zakharov spectrum. 



1.5 Main results and predictions 



1.5.1 Alfven wave turbulence 

We start to summarize the results of wave turbulence in the incompressible MHD 
case for which the wave kinetic equations are singular in the sense that the principal 
value terms remain except for the Elsasser energies E ± . As it was explained, the ori- 
gin of this particularity is the pseudo-dispersive nature of Alfven waves which are a 
uniqu e case where wave turbulence theory applies (see also IZakharov and Sagdeevl . 
1970 ). In the limit of strongly anisotropic turbulence for which k± ^> fcii, the wave 
kinetic equations (|1.38p simplify. Then, the Elsasser energy spectra of the trans- 
verse fluctuations (shear-Alfven waves) satisfy the asymptotic integro-differential 
equations 

dE ± (k±,k\ ] ) _ 



dt 



(1.40) 



7T£ 

"Bo" 



>sin0 -^£T(q ± ,0) [k ± E ± (p ± ,k ll ) - p {k ^ k\\)] dp ± dq ± . 
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where A defines the domain of integration on which kj^ = px + q_i_ , <t> is the angle 
between kj^ and pj_ and the angle between k^ and qj_. A fundamental property 
of Alfven wave turbulence appears on equations (| 1.400 : the evolution of the spectra 
E^ is always mediated by interaction with the slow mode qu = 0. This is a direct 
consequence of the resonance conditions discussed in Section 11.3.2.31 Then, the 
parallel wavenumber dependence does not affect the nonlinear dynamics: in other 
words, there is no transfer along Bo and fc|| can be treated as an external parameter. 
Therefore, we introduce the notation E/(k±,kn) = E/(k±)f±(ku) where f± are 
arbitrary functions of fc|| given by the initial conditions and with the assumption 
/±(0) = 1 we obtain 

dE ± (k±) 



dt 

■XS 2 f f n , . „ k 



(1.41) 



cos 2 0sin0— E*(q ± ) [k x E ± (p ± ) - p±E ± (k 1 _)] d P± dq ± , 
-do JJa Qx 

which describes the transverse dynamics. The exact finite flux solutions of equations 

(|1.41[) for the stationary energy spectra are 

E ± (k 1 _)^k n A ± (1.42) 

with 

n++n- = -4. (1.43) 

Additionally, the power law indices must satisfy the condition of locality (i.e. the 
condition on the power law solutions for which the wave kinetic equations remain 
finite) 

-3<n±<-l. (1.44) 

It is important to recall that wave turbulence is an asymptotic theory which must 
satisfy a condition of application. In our case the transfer time must be significantly 
larger than the Alfven time which means in terms of wavenumbers that k\\ ^> e 2 k± 
(and that at small transverse scales turbulence becomes strong). At this point 
a comment has to made about the Kolmogorov-Zakharov solutions (11.430 . These 
solutions imply the contribution of the wavevector q which is by nature a slow mode 
(<7|l = 0) whereas the other contributions k and p imply wave modes (k» — p|| > 0). 
In Fig. 11.31 it is assumed that E + = E~ and that the dynamics of the slow mode 
is given by the wave mode hence the stationary solution n + = n_ = 2. However, it 
may happen that the slow mode has its own dynamics which can belong to strong 
turbulence. Since relation (|1.43|) is exact, if the wave mode energy has a —5/3 scaling 
then the wave turbulence s pectrum should be —7/3. Although this discussion was 



given in the original paper (jGaltier et all l2000[) it was not considered seriously until 



the m ost recent direct numerical simulations ( Bigot et al. , 2008tJ : Bigot and Galtier 



20111 ) reveal that this situation may happen. We will come back to this point in 



the conclusion. 
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The locality of interactions (in Fourier space) is an important issue in MHD 
turbulence. According to some recent works in isotropic turbulence, nonlinear in- 
teractions seem to be more non local in MHD than in a pure hydrodynamics in the 
sense that the transfer of energy from the vel ocity field to the mag netic field may 



be a highly non local process in Fourier space (jAlexakis et all 120051 ) . The situation 



ragr 



is different when we deal with anisotropic turbulenc e: in this case interac tions (be 



tween perpendicular wavevectors) are mainly local ([Alexakis et all 120071) . In wave 



turbulence, the condition (|1.44p has to be satisfied to validate the exact power law 
solutions and avoid any divergence of integrals in the wave kinetic equations due 
to nonlocal contributions. In practice, numerical simulations of th e wave kinetic 



equat ions have clearly shown that the solutions (|1.44[) are attractive (jGaltier et al 
2000l) . 



Recently, it was realized that the wave kinetic equations found in the anisotropic 
limit may be recovered without the wave number condition k± ^> k\\ if initially only 
the shear- Alfven waves were considered (jGaltier and Chandranl . 120061 ). Therefore, 
the finite flux solution may be extended to the entire wavenumber space which 
renders its de tection easier. This idea w as tested successfully against numerical 
simulations by Perez and Boldvrev (2008) who found a spectral signatures of Alfven 
wave turbulence. 



1.5.2 Compressible MHD 

We turn now to the compressible regime for which two limits have been analyzed. 
The first case (case I) is the one for which the plasma (thermal) pressure is assumed 
to be small as compared to the magnetic pressure (small j3 limit) and where three- 
wave interaction processes that do not involve slow waves are neglected ([Kuznetsovl . 
20011 ). In the second case (case II) for which we still have <C 1, slow waves are 
neglected and a constant density is imposed (jChandranl . 120051 ) . In both situations 
the general finite flux solutions are not obvious to express. 

In case I, when only interactions between Alfven and slow waves are kept, a wave 
energy spectrum in ~ k^k^ is found which corresponds to a (finite) constant 
energy flux solution. It is claimed that the addition of the interactions with the fast 
waves will lead to the same solution since the dynamics tends to produce strongly 
anisotropic distributions of the waves concentrated in k-space within a narrow-angle 
cone in the Bo direction (k± <C fcii). Under these conditions, the fast waves coincide 
with the Alfven waves. 

In case II, the general wave kinetic equations do not allow us to find exact power 
law solutions. However, when only Fast-Fast-Fast interactions are kept it is possible 
to find a finite flux solution for the fast wave ID energy spectrum which scales as 
~ /(0)fc" 3/2 , where f{6) is an arbitrary function of the angle 9 between the wave 
vector k and the uniform magnetic field Bo- Numerical simulations of the general 
wave kinetic equations are made to find the behavior according to the angle 0. A 
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solution close to fcj 2 is found for the Alfven wave 2D energy spectrum (for different 
fixed fe|i) when k± ^> k\\. The fc-spectra plotted at 6 = 45° reveal a fast wave 
modal energy spectrum in ~ fc~ 3 / 2 and an steeper Alfven wave spectrum, while 
for a s mall angle (0 = 7 .1°) both spectra follow the same scaling law steeper than 
k' 3 / 2 dChandranlHooil) . 



1.5.3 Whistler wave turbulence 

The electron MHD equations in the presence of a strong uniform magnetic field 
Bo exhibit dispersive whistle r waves. The wave turbulence regime was analyzed 
in the incompressible case by Galtier and Bhattachariee ( 20031 2005 ) who derived 
the wave kinetic equations by using a complex helicity decomposition. The strong 
anisotropic (kj_ k») finite flux solutions correspond to 



E(k x M) 



,-5/2,-1/2 
Hj_ rvii 



for the magnetic energy spectrum, and 

ff(fcx,fc|i)~fcl r/2 fc|7 1/2 



(1.45) 



(1.46) 



for the magnetic helicity spectrum. As for the other cases presented above, a direct 
cascade was found for the energy. In particular it was shown that contrary to 
MHD, the wave kinetic equations which involve three- wave interaction processes 
are characterized by a nonlinear transfer that decreases linearly with k\\ ; for fc|| = 0, 
the transfer is exactly null. Thus the 2D modes (or slow modes) decouple from 
the three-dimensional whistler wa ves. Su c h a decoupling is found in a var iety of 
problems like rotating turbulence ( Galtier . 2003 ; Sagaut and Cambonl . 20081 ). 



1.5.4 Hall MHD 

The last example exposed in this chapter is the Hall MHD case which incorporates 
both the standard MHD and electron MHD limits. This system is much heavier to 
analyze in the r egime o f wave t urbulence and it is only recently that a theory has 
been proposed (Galtier, 2006bf) . The general theory emphasizes the fact that the 
large scale limit of standard MHD becomes singular with the apparition of a new 
type of terms, the principal value terms. Of course, the large scale and small scale 
limits tend to the appropriate theories (MHD and electron MHD wave turbulence), 
but in addition it is possible to describe the connection between them at intermedi- 
ate scales (scales of the order of the ion inertial length d,). For example, moderate 
anisotropy is predicted at these intermediate scales whereas it is much stronger for 
other scales. It is also interesting to note that the small scale limit gives a system of 
equat ions richer than the pure electron MHD system ([Galtier and Bhattacharjed . 
20031 ) with the possibility to describe the sub-dominant kinetic energy dynamics. An 
exact power law solution for the kinetic energy spectrum is found for ion cyclotron 
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wave turbulence which scales as 

E(k±,k {] ) ~ fc; 5/2 fc[ 1/2 . (1.47) 

To date, the wave kinetic equations of Hall MHD have not been simulated numer- 
ically even in their simplified form (when helicity terms are discarded). It is an 
essential step to understand much better the dynamics at intermediate scales. 



1.6 Conclusion and perspectives 
1.6.1 Observations 

Waves and turbulence are two fundamental ingredients of many space plasmas. The 
most spectacular illustration of such characteristics is probably given by the obser- 
vations of the Sun's atmosphere with the orbital solar observatory Hinode/JAXA 
and more recently by SDO/NASA launched in 2010. For the first time, detection of 
Alfven waves is made through the small oscillations of many thin structures called 
threads (see Fig. I1.4|) . In the meantime the highly dynamical nature of coronal loops 
is revealed by non-thermal velocities detected with spectrometers. These findings 
are considered as a remarkable step in our understanding of the solar coronal dy- 
namics. Nowadays, it is believed that Alfven w aves turbulence is a promising model 



to understand the heating of the solar corona (jvan Ballegooiien et all 120111 ). 

The interplanetary medium is another example of magnetized plasma where 
waves and turbulence are detected. In this framework, the origin of the so-called 
" dissipative range" , i. e. the extension of the turbulent inertial range beyond a frac- 
tion of hertz, is currently one of the main issue discussed in the community. Al- 
though a final answer is not given yet, wave turbulence is a promising regime to 
understand the inner solar wind dynamics in the sense that it gives exact results 
in regards to the possible multiscale behavior of magnetized plasmas as well as the 
intensity of the anisotropic transfer between modes. 

The main feature of magnetized plasmas in the regime of wave turbulence is 
the omnipresence of anisotropy and the possibility to have different spectral scaling 
laws according to the space direction. To achieve a proper comparison between 
observational data and theoretical predictions, not only in situ measurements are 
necessary, but multipoints data have also to be accessible. It is at this price that 
the true nature of magnetized plasmas will be revealed. The magnetosphere is the 
first medium whe re it is possible to per form such a comparison thanks to the Clus- 



ter/ESA mission (jSahraoui et al.l . 120061 ) . In the case of the Jupiter's magnetosphere 



the large scale magnetic field is relatively strong. Actually, the firs t indirect evi- 



dence of Alfven wave turbulence was reported in ISaur et al.l (|2002f ) by using five 
years set of Galileo/NASA spacecraft magnetic field data. (Since the work was 
based on a one point analysis a model was used to differentiate the perpendicular 
and parallel spectra.) 



January 9, 2012 1:7 



World Scientific Book - 9.75in x 6.5in 



Galtier'book 



Wave turbulence in magnetohydrodynamics 



27 




Fig. 1.4 Schematic view of Alfven wave turbulence on the Sun: coronal loops act as a reso- 
nant cavity for Alfven wave packets which are transmitted from the deeper layers of the solar 
atmosphere. 



Perhaps the first direct evidence of Alfven wave turbulence comes from solar 
(photospheri c) magnetic field o bservations above active regions. Using data from 
SOHO/ESA. lAbramenkol (|2005l ) observed that the magnetic energy spectra reduced 
to its radial component (along Bo) have a transverse scaling significantly different 
from the Kolmogorov one with power laws generally between —2 and —2.3. Since 
the magnetic field can only be detected in a thin layer (the temperature is too 
high in the chromosphere and in the corona to measure the magnetic field with 
the Zeeman effect), the magnetic spectra can be interpreted as a measure at a 
given fen > (in practice a small average is made in feii since the data comes from 
slightly different altitudes). As we know from the Alfven wave turbulence theory 
the parallel fluctuations follow the same dynamics as the perpendicular (shear- 
Alfvcn) fluctuations and the addition of the power law indices of the slow and wave 
mode spectra is —4. Then, these unexpected spectra may be seen as the direct 
manifestation of Alfven wave turbulence if the slow mode contribution scales like 
Kolmogorov (see the discussion on Section [T33J- If this interpretation is correct it 
is probably the most important achievement of the Alfven wave turbulence theory. 



1.6.2 Simulations 

Numerical simulation is currently the main tool to improve our knowledge on wave 
turbulence in magnetized plasmas since we are still limited by the observational 
(single point) data. Two types of simulations are available: the simulation of the 
wave kinetic equations and the direct numerical simulation of the original MHD- 
type equations. In the former case, it is a way to find for instance the spectral 
scaling laws when the wave kinetic equations are too complex to provide us the 
exact solutions after application of the usual conformal transform. An example is 
given by compressible MHD: the numerical simulations revealed a variation of the 
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power law energy spectrum with the angle between k and Bo ( Chandran , 20051) . 
Simulations may also be useful to investigate wave turbulence when an ex ternal 
forcing is applied like in incompressible MHD (jGaltier and Nazarenkol 20081 ). 

For direct numerical simulation, the challenge is slightly different: indeed, in 
this case the main goal is the measure of the transition between strong and wave 
turbulence, and thus between isotropic and anisotropic turbulence. The former 
regime has been extensively studied since more than three decades whereas the 
latter is still a young subject. The main topic of such a simulation is also to find 
general properties that could help us to understand the single point measurements 
made in natural plasmas. We arrive here at the heart of current issues in wave 
turbulence. One of the most important points emphasized by recent direct nu- 
merical sim ulations in incompr essible MHD is the coexistence of wave and strong 



turbulence ( Bigot et al. . 2008bl) . This characteristic should not be a surprise since 



basically wave turbulence is a perturbative theory which must satisfy conditions 
of applicability. In this case, the slow mode (<7|| = 0) may evolve differently from 
the wave modes (fc|hP|i > 0) since the former case may be characterized by strong 
turbulence and the latter by wave turbulence. This distinction is fundamental and 
not really taken into account in the community: Alfven wave turbulence will not 
be fully revealed in space as well as simulated plasmas as long as strong and wave 
turbulen ce will not be separated. This point is illustra ted with direct numerical sim- 
ulations (| Bigot et all l2008bt iBigot and Galtierl . 120111 ) where steeper energy spectra 
may be found for the wave modes (between —2 and —7/3) compared to the slow 
mode (down to —5/3). The result seems to depend on the relative intensity of B : 
a strong magnetic background tends to increase the energy of the slow mode which, 
then, may evolve independently (the critical value of Bq seems to be around 10 
times the r.m.s. fluctuations). It is also found that the wave turbulence spectra are 
generally visible at a fixed fcii > and may be hidden in a —5/3 scaling if a summa- 
tion over fcii is performed (this result could be due to the limited space resolution). 
The same investigation clearly shows that - as expected- an equipartition between 
the kinetic and magnetic energies is reached in the wave modes whereas the mag- 
netic energy is dominant in the slow mode. Note th at this Alfven wave turbule nce 
regime was also numerically found independently bv lPerez and Boldvrevl (|2008l ). 



1.6.3 Open questions 

In the light of such recent results the following new questions may be addressed: 

• Docs the —7/3 spectrum more universal than —2 in Alfven wave turbulence ? 

• Is the slow mode (qu = 0) intermittent ? and then does it produce intermit- 
tency in the slave wave modes (fc||,P|| > 0) ? 

• Can we find direct signatures of Alfven wave turbulence in solar magnetic 
loops with three-dimensional magnetic field data ? 
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• Does Alfven wave turbulence the universal regime for stellar magnetic loops ? 

• Does the —5/3 energy spectrum observed in the solar wind correspond to a 
bias ? and then what is the wave vector energy spectrum ? 

• Is the small energy ratio E u /E b ~ 0.5 found in the local solar wind 
(jBruno and Carbond . 120051) due essentially to the slow mode contribution in which 
case the wave modes may follow equipartition ? 

• Are the absence of intermitten cy and the —2.5 m agnetic spectrum observed 
in the solar wind dispersive regime (jKivani et all 120091 ) the signatures of whistler 
wave turbulence ? 



Future missions like Solar Orbiter from the European Space Agency (scheduled 
for 2017) will certainly help us to answer these questions and make significant 
advances on plasma turbulence. 
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